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Elastic pd scattering is studied within the Glauber-Sitenko theory based on the single- and double 
pN scattering mechanisms. The full spin dependence of the elementary pN scattering amplitudes 
is taken into account and both the S- and D-wave components of the deuteron are considered. The 
treatment of the spin dependence is done in a (properly modified) formalism developed recently by 
Platonova and Kukulin for the pd — > pd scattering process. Predictions for differential cross sections 
and the spin observables Ay, A%, A xx , A yy are presented for antiproton beam energies between 50 
and 300 MeV, using amplitudes generated from the NN interaction model developed by the Jiilich 
group. Total polarized cross sections are calculated utilizing the optical theorem. The efficiency of 
the polarization buildup for antiprotons in a storage ring is investigated. 
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I. INTRODUCTION 

The present investigation is motivated by the plans of 
the PAX collaboration [l[ to measure the transversity of 
the proton (antiproton) in double-polarized Drell-Yang 
processes at an upgrade of the FAIR facility in Darm- 
stadt. In order to achieve this aim an intense polarized 
beam of antiprotons is required. A possibility to over- 
come the experimental challenge to obtain such a po- 
larized beam is seen in elastic scattering of antiprotons 
off a polarized 1 H target in rings [2]. Analogous exper- 
iments performed for the proton case by the FILTEX 
collaboration [3[ at 23 MeV and a recent COSY study 
where protons were scattering off a polarized hydrogen 
at 49 MeV 0| showed that indeed a polarized (proton) 
beam can be achieved via the so-called spin-filtering ef- 
fect. According to theoretical interpretations |5j of the 
data [3, |j| , the polarization buildup effect appears solely 
due to the hadronic interaction of the incoming proton 
with the target. 

Whereas the spin dependence of the nucleon-nucleon 
(NN) interaction is very well known at the considered 
energies, that allows one to calculate reliably the spin- 
filtering effect for protons, there is practically no corre- 
sponding information for the antiproton-nucleon (NN) 
interaction. For this reason a test experiment for the 
spin-filtering effect in the antiproton-hydrogen interac- 
tion is planned at the AD ring at the CERN facility [f|Q- 

In view of the unknown spin dependence of the pN 
interaction, the interaction of antiprotons with a po- 
larized deuteron is also of interest for the issue of the 
antiproton polarization buildup. This option was dis- 



cussed in a previous paper by us [S]. In that work the 
single-scattering approximation was used for the calcu- 
lation of the polarized total pd cross sections for ener- 
gies in the region 50-300 MeV. The spin dependence of 
the elementary pN amplitudes was taken into account in 
this approximation only in collinear kinematics using the 
NN interaction model of the Jiilich group [9T4T21]. The 
pN double-scattering effects were only accounted for in 
the computation of the unpolarized total and differential 
cross sections and found to be in the order of 10-15% 
8]. Spin observables for pd elastic scattering and shad- 
owing effects (double scattering) in polarized total cross 
sections were not considered in that work. A calculation 
of such observables, including double-scattering effects, is 
the main aim of the present work. Indeed corresponding 
results are certainly desirable, specifically in view of the 
prospect that a discrimination of existing models of the 
pN interaction could be feasible on the basis of a com- 
parison with expected data Q- Antiproton polarization 
buildup in the context of pd scattering was also studied 
in a recent work by Salnikov (T3 |. utilizing amplitudes 
taken from the Nijmegen pp partial wave analysis fl4| . 

In the present paper we consider elastic pd scatter- 
ing within the Glauber-Sitenko theory of multi-step scat- 
tering [l5|], taking into account the full spin-dependence 
of the elementary pN scattering amplitudes. There are 
several studies of the accuracy of the Glauber-Sitenko 
theory in the literature [Rj-[2l| which demonstrate that 
corrections to the eikonal approximation, which is the 
basis of this theory, are small in the region of interme- 
diate energies about ~ 1 GeV. The reliability of the 
Glauber-Sitenko approach at intermediate energies was 
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studied recently in Ref. [22j via a comparison with rigor- 
ous Faddeev calculations for the case of identical spinless 
bosons interacting by means of a simple Malfliet-Tjon 
interaction potential. The results of Rcf. [22] for such a 
"bosonic" nd system show that even at rather low ener- 
gies, 100 ~ 200 MeV, the difference between the Faddeev 
and the Glauber-Sitenko calculations is just about 10- 
15% for the total cross section. Rather good agreement 
was found also for the differential cross section in the for- 
ward hemisphere, excluding the region of the diffraction 
minimum [221 ] . 

For antiproton-nucleus scattering the Glauber-Sitenko 
theory can be applied at lower energies [H, [24| as com- 
pared to the proton-nucleus reaction. The amplitude 
for elastic scattering of antiprotons off nuclei is strongly 
peaked in forward direction due to strong annihilation ef- 
fects in the pN interaction, supporting the applicability 
of the eikonal approximation (see Ref. Q and reference 
therein) . In the present work the S- and D-wave compo- 
nents of the deuteron and both the single- and double pN 
scattering mechanisms are taken into account. The treat- 
ment of the spin dependence is based on a proper modifi- 
cation of the formalism developed recently by Platonova 
and Kukulin [25|, [26[ for pd elastic scattering. In their 
papers the formalism was succesfully applied for describ- 
ing spin observables of the pd — > pd process at 250-1000 
MeV. An independent confirmation of the findings of [25| 
was reported recently by us 28] . 

The spin dependence of the pd — ► pd amplitude is very 
similar to that for pd — > pd scattering, except for the 
contribution of the charge-exchange channel pp O fin 
which, however, can be taken into account straightfor- 
wardly. In the present work we consider the differential 
cross section and the spin observables Ay, AZ, A xx , A yy 
for the pd — > pd process. Those observables are evaluated 
here for antiproton beam energies from 50 to 300 MeV 
employing pN amplitudes generated from the Jiilich NN 
model ll]. The total polarized cross sections o\, a%, and 
(73 are calculated on the basis of the optical theorem. We 
also investigate the efficiency of the polarization buildup 
for antiprotons in a storage ring. Here Coulomb effects 
are taken into account within the formalism described in 
Ref. @. 

The paper is structured in the following way: In Sect. 
II we introduce briefly the used formalism. In particular, 
we point out the differences that occur between its ap- 
plication to the pd- and to the pd systems. In Sect. Ill 
results for pd scattering are presented. First we discuss 
the issue of the applicability of the Glauber-Sitenko the- 
ory. Specifically, we assess the angular range for which 
differential observables can be reliably calculated within 
this approach. Then predictions for the differential cross 
section and the spin observables Ay, A?, A xx , A yy are 
given based on elementary pN amplitudes taken from the 
Jiilich NN models A and D. Finally, our results for the 
total pd cross section are provided, including those for 
the polarized case. We also provide predictions for the 
polarization efficiency which is the decisive quantity for 



the spin-filtering method. The paper ends with a short 
Summary. Relations between amplitudes and considered 
observables are given in an Appendix. 



II. ELEMENTS OF THE FORMALISM FOR 
pd -> pd SCATTERING 

For the pd —> pd process we use the formalism devel- 
oped in Ref. @ for the process pd — > pd, taking into 
account the specific differences that arise for the pd col- 
lision. Within the Glauber-Sitenko theory [l5| the scat- 
tering matrix for elastic pd scattering is given by the 
following matrix element 



M P M) = (/|M(q,s)|t> 



(1) 



evaluated between definite initial |i) and final |/) states of 
the pd system. Here q is the momentum transferred from 
the initial to the final deuteron in the pd collision and s 
is the impact parameter of the projectile. In Eq. ([T]) and 
in the following we suppress the dependence of Mpd on 
the total energy to simplify the notation. The transition 
operator M can be written as 

M(q,s) = 

exp (-iq ■ s)Alpp(q) + exp (- -iq ■ s)M pn (q) + 



27r 3/2 j exp (* q/ ■ s ) [ M pp(qi) M p«(q2) + 

Mp„(qi)Mp P (q 2 ) - Mp P ^ fi „(qi)M ftn ^p P (q 2 ) d 2 q'. 

(2) 

Here Mp_/v(q) (N — p,n) is the pN scattering matrix, 
->pp) is the scattering matrix of the charge- 
exchange process pp — > fin (nn — > pp), and qi = q/2 — q' 
and q 2 = q/2 + q' are the transferred momenta in the 
first and second pN collision, respectively, in the double- 
scattering terms. The first two plane- wave terms corre- 
spond to the single-scattering mechanism, while the three 
terms in the integral represent the double-scattering 
mechanism. These transitions correspond to the dia- 
grams depicted in Fig. [TJ It is assumed in Eq. ^ that the 
deuteron wave function in the matrix element of this op- 
erator does not contain the isospin part explicitly. There- 
fore the last term in Eq. ([2]) appears with a negative sign 
coming from the product of the isospin factors at the dpn 
vertices. 

The elementary scattering matrix for elastic pN scat- 
tering is given by the following expression: 

M pN = A N + [C n cti + C' N (T2) ■ n 

+ B N {(T 1 -k)((T2-'k) 

+ (G N - H N )(ax ■ n)(<7 2 • n) 

+ (G JV +-ffjv)(o-i-q)(cr 2 -q) . (3) 

Here tj\ (0*2) is the Pauli matrix acting on the spin of 
the p (N) states, N = p,n. The unit vectors are defined 
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by k = (ki + k/Vlki + k/l, q = (k* - k/)/|ki - k/|, and 
n = [k x q], where k^ (k/) denotes the momentum of 
the incident (outgoing) antiproton. The one for charge- 
exchange, Alpp^nn, has the same spin structure as given 
in Eq. ([3|). We denote the corresponding amplitudes on 
the right-hand side of Eq. ([3]) by A c , C c , etc. in the 
following for simplicity reasons. 

The amplitudes in Eqs. ([IJ-© are normalized as in 
Ref. [25[ and related to the corresponding differential 
cross sections via 

J = \s P M pd M+ d , (4) 

^ = \s P M pN M+ N , (5) 

where t = — q 2 is the squared four-momentum transfer in 
the pd or pN systems, respectively. The factor in front 
of the integral in Eq. © differs from that in the original 
Glauber theory [1 51 ] due to a different normalization of 
the invariant scattering matrices M used in Eqs. (HJ-©. 

In general there are 36 transitions in the pd — » pd 
process for the different spin states. When accounting 
for rotational invariance and invariance under parity- 
and time-reversal transformations, the number of tran- 
sition matrix elements for the pd — > pd process re- 
duces to twelve independent complex amplitudes A4 (i = 
1, . . . , 12), which can be introduced in the same way as 
for the pd —> pd process [1,[25[ (see also the Appendix). 
The amplitudes A4 can be written in the form 

Mi) = 4 s) (q) 

+ 2^75 / d\'(AfM)- 4 c) (q,q')).(6) 

Here we introduced the scalar amplitudes for single scat- 
tering (Fig. [Th), A\ s \q), double scattering (Fig. [lb), 
A[ d \q,q'), and double-scattering involving charge ex- 
change (Fig.QJ), (q, q') using the same notations as 
in Ref. [25|. 

The amplitudes A^f' 11 '^ can be expressed via the elastic 
form factors of the deuteron corresponding to the transi- 
tions S^S,S^D,D^S,D^D, from the initial to 
the final deuteron state, and the elementary amplitudes 
of pN scattering. The explicit expressions can be found 
in Ref. [25| in the appendices A and B and, therefore, 
we refrain from reproducing those lengthy formulae here. 
Specifically, the formulae for the single-scattering mecha- 
nism Af^ (Fig. [IJi) are the same as those for the pd — » pd 
processs given in Ref. [25[ in Table I with the replacement 
of the elementary amplitudes of pN elastic scattering by 
the corresponding ones for pN . However, the integral in 
Eq. ^ , accounting for the double-scattering mechanisms 
differs from that in Eq. (A.l) in Ref. [25[ by two aspects. 
First, the term A^ in Ref. [25| has an additional fac- 
tor 2. Second, the permutations n O p are taken into 
account here for A\ s \q) and A[ d \q, q'), as in Ref. [25$ . 
but they do not apply for A[ c \q, q'). 




d d 
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FIG. 1: Mechanisms for pd elastic scattering: single scattering 
(a), double scattering (b), and charge-exchange (c). 

In order to explain the reason for these two differences 
between the pd — > pd and pd — > pd amplitudes, let us 
discuss first the pd — > pd process following Ref. (25J. In 
the pd — > pd process there are three double-scattering 
amplitudes, which we denote symbolically as 

M pp (qi) x M pn (q 2 ) + M p „(qi) x Af pp (q 2 ) 

(qi) x Af rlp ^p„(q 2 ) . (7) 

The small-angle charge-exchange scattering matrix 
M pn ^ np can be rewritten in terms of the small-angle 
scattering matrices M pp ^ pp = M pp and M pn ^ pn = M pn 
in the form M pn ^ np (q) = M np ^ pn (q) = M pp (q) - 
M p „(q). (The double-charge-exchange scattering matrix 
vanishes in the approximation M pp (q) — M pn (q), i.e. if 
one disregards the isospin dependence of the pN ampli- 
tude.) Therefore, the third term in Eq. ((7]) can be written 
as 

(qi) x M„ p ^ pn (q 2 ) = 
Afpp(qi) x Mpp(q 2 ) + Mp„( qi ) x M p „(q 2 ) 
-Mpp( qi ) x A/p„(q 2 ) - Afpn(qi) X M p p(q 2 ) . (8) 

The term M pn ^ np (c\x) x M np ^ pn {c\2) enters the full pd 
amplitude with opposite sign with respect to the first two 
terms given in Eq. (JJJ due to the additional permutation 
pf>nin the deuteron vertex. Furthermore, the first and 
the second terms on the right hand side of Eq. ((SJ con- 
stitute the term called in Ref. (2f| the " charge-exchange" 

(c) 

amplitude and denoted by A\ . The third and the fourth 
terms from Eq. (|8]) can be absorbed into the first and sec- 
ond terms in Eq. (J7J, respectively, and this leads to the 

factor 2 in front of the A\ d ' amplitude in [25[ . 
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Coming back to the process pd — ► pd, one should 
mention that in this case there are also three double- 
scattering terms in Eq. ([2]), namely 

Mppfax) x M Pn {q2), Mp„(qi) x Afp P (q 2 ), 

Mp P _j. Sn (qi) x M fin _>p P (q 2 ) . (9) 

The first two terms constitue the amplitude in 

Eq. ©. The detailed formulae for A { ? ] (i =_L . . . , 12) 
coincide with those given in Table II of Ref. |25[ for the 
pd — » p<i, where each term (An, Bn, Cat, CJy, Gat, Hn) 
from the piV — > pN amplitude (N — p, n) has to be re- 
placed by the corresponding term of the scattering matrix 
MpN^pM given in Eq. ([3]). 

The last term in Eq. ^ gives rise to the A| c ' ampli- 
tude in Eq. and, like for the pd — > pd process, en- 
ters the full amplitude also with the opposite sign with 
respect to the first two amplitudes in Eq. @ for the 
same reason as in the pd — > pd. In the present cal- 
culation the term with the charge- exchange amplitudes 
pp o fin is not expressed through the small-angle scat- 
tering amplitude, as in the pd — ¥ pd case, but calcu- 
lated straightforwardly and, therefore, does not lead to 
the factor 2 in front of in Eq. ([6]) in contrast to 

Ref. H3- The formulae for A[ c} (i = 1, . . . , 12) are the 
same as those given in Table III of Ref. [25[ , where each 
product A„(q 2 )A„(qi), etc., has to be replaced by the 
corresponding terms A c (q 2 )A c (qi), etc., of the charge- 
exchange scattering matrix Mp pi ^n n . Finally, we want 
to emphasize that contrary to Ref. [25| the total double- 
scattering term enters Eq. © with positive sign. The 
relations of these amplitudes Ai (q) with the spin observ- 
ables considered in the present paper are given in the 
Appendix. 

III. NUMERICAL RESULTS AND DISCUSSION 
A. General remarks and differential observables 

Earlier studies [H, HH HH and also our previous cal- 
culations [1, H3] were all done within the spinless ap- 
proximation for the elementary pN amplitude Mpjv, i.e. 
using only An from Eq. ([3]) parametrized in the Gaus- 
sian form, and restricted to using only the 5*- wave part of 
the wave function of the target nucleus. Those investiga- 
tions suggested that the Glauber-Sitenko theory allowed 
one to explain the differential cross sections in the for- 
ward hemisphere and the total unpolarized cross section 
for the reactions pd, p 3 He and p 4 He even at rather low 
energies, i.e. down to 20-50 MeV of the incident an- 
tiproton. For p 4 He elastic scattering the first and the 
second diffraction peaks were explained by these calcula- 
tions [13] ■ On the other hand, attempts to describe the 
second peak observed in pd elastic scattering at 179 MeV 
[3l( | were not quite that successful as documented in some 
papers [U, Elf- 



in the present calculation we include the S- and D- 
wave components of the deuteron wave function and 
we keep the full spin-dependence of the pN amplitude 
as given in Eq. ([3]). For the former we use the wave 
function of the CD Bonn potential as parameterized in 
Ref. [32| • For generating the latter we employ two mod- 
els developed by the Jiilich group, namely A(BOX) in- 
troduced in Ref. @ and D described in Ref. |Tl]. Re- 
sults for the total and integrated elastic (pp) and charge- 
exchange (pp — ¥ fin) cross sections and also for angular 
dependent observables for both models can be found in 
Refs. H,[ll| and, specifically, spin-dependent observables 
in [liq . Model A as well as D provide a very good over- 
all reproduction of the low- and intermediate energy NN 
data as documented in those works. 

An exemplary result demonstrating the role of the 
single-scattering (SS) and double-scattering (DS) mecha- 
nisms is presented in Fig. [5] There, cross sections based 
on the SS- and the DS mechanisms and their coherent 
sum (SS+DS) are displayed separately. One can see that 
the SS mechanism alone fails to explain the forward peak. 
However, the sum SS+DS describes it rather well. Ob- 
viously, the DS mechanism, neglected in Ref. Q in the 
calculation of the spin-dependent total cross sections, has 
a sizable influence even in the region of the forward peak. 
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FIG. 2: Differential cross section of elastic pd scattering at 
179 MeV versus the transferred momentum squared. Re- 
sults based on single-scattering (dash-dotted line), double- 
scattering (dashed) and the full (solid) Glauber-Sitenko mech- 
anisms are shown. For the calculation the pN amplitude pre- 
dicted by the Jiilich model D is utilized. The data points are 
from Ref. [3J. 

When accounting for the spin dependence of the pN 
interaction one has to address also the reliability of the 
employed Glauber-Sitenko approach, specifically, as far 
as the angular range is concerned. In contrast to the 
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FIG. 3: Differential elastic pp (a) and pn (b) cross sections 
at 179 MeV as predicted by model D (solid lines). The 
dashed line shows the result with the smooth cutoff with 
go = 0.21 GeV/c and v = 10, corresponding to 6 p c N m . = 105°, 
as explained in the text. 



spin-independent part of the amplitude Eq. given by 
Am (N — p,n), most of the other amplitudes that give 
rise to the spin dependence (-Bat, Cat, C' N , Gn, Hn in 
Eq. ©) do not exhibit a well-pronounced diffractive be- 
haviour for antiproton beam energies 50-200 MeV, i.e. 
they do not decrease rapidly with increasing center-of- 
mass (cm.) scattering angle 9 c . m .. In fact, some of these 
amplitudes even increase with increasing 9 c , m . and their 
magnitude is larger at 8 c . m . > 90° than at 9 c . m . < 90°. 
The typical behaviour of the differential pN cross sec- 
tion can be seen in Fig. [3] for 179 MeV. This behaviour 
is basically the same for both Julich models A and D 
and, therefore, we consider only model D in the follow- 
ing. The visible rise of the cross section for backward 
angles is partly due to those spin-dependent pieces of 
the pN amplitude and is reflected in the corresponding 
pd results by the appearance of a second, broad peak at 
large transferred momenta q 2 ~ 0.35 (GeV/c) 2 coming 
from the SS contribution alone, as is seen in Fig. [5J As 
a consequence, the second ("diffraction") peak that ap- 
pears in the full calculation, including now the SS and 
DS mechanisms, originates not only from the interfer- 
ence between the SS and DS amplitudes, as is usually 
the case for typical diffractive scattering, but is also re- 
lated to that backward peak seen in the cross section 
of the (elementary) pN reaction. We confirmed in test 
calculations within the SS mechanism, utilizing just the 
amplitudes An in Eq. © or a standard Gaussian- type 
representation of the total pN amplitude in the spinless 
approximation, that then the differential cross sections 
do not demonstrate this behavior. 

One should note, that the Glauber-Sitenko diffraction 
theory of multi-step scattering is not suitable for tak- 
ing into account backward scattering in the elementary 
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FIG. 4: Differential cross section of elastic pd scattering at 179 
MeV versus the cm. scattering angle. Results of our full cal- 
culation (including the SS+DS mechanisms) are shown based 
on the NN models model A (green/grey) and D (red/black). 
The bands represent the sensitivy to variations of the large- 
angle tail of the pN amplitudes as discussed in the text. The 
data points are taken from Ref. [311 ]. 
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FIG. 5: Differential cross section of elastic pd scattering at 
282 MeV versus the cm. scattering angle. Same description 
of curves as in Fig. [J] 

hadron-nucleon collision, because its basis is the eikonal 
approximation. The deuteron elastic form factor sup- 
presses the contribution of the pN amplitudes at large 
angles, but at the considered low and moderate energies 
this effect is not as strong as for higher energies (several 
hundred MeV) for the same scattering angle. 

The authors of Refs. [25|, [26[ approximate all ampli- 
tudes in Eq. © as a sum of Gaussians, e.g., 

A N(q) = ^c AN . J e^p{-d AN jq 2 ), etc., (10) 
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FIG. 6: Differential cross section of elastic pd scattering at 50 
MeV versus the cm. scattering angle. Same description of 
curves as in Fig. 2] 



where ca n j and d,A N ,j are parameters fixed by a fit to 
the original NN amplitudes. With such a series of Gaus- 
sians they are able to reproduce the original amplitudes 
quite well in the forward hemisphere, i.e. up to angles of 
8 c .m. = 80° — 90°, as reported in [2g|. For larger angles 
the differential NN cross section based on those Gaus- 
sians drop off very quickly, in accordings with the re- 
quirements by the Glauber-Sitenko model, as shown in 
Fig. 4 of Ref. (26|. However, it remains unclear how and, 
specifically, how fast the individual (parameterized) am- 
plitudes drop off. In particular, issues like the influence 
of the backward tail on the results or the related question 
of up to which angles one can trust the predictions for 
pd observables are not addressed in Ref. |25j . 

In the following we investigate this subject in the con- 
text of elastic pd scattering. For this purpose we perform 
various calculations of this reaction within the Glauber- 
Sitenko theory in order to scrutinize the sensitivity of the 
considered pd observables to the large- angle region of the 
employed elementary pN amplitudes. Clearly, any such 
sensitivity is in contradiction with the assumptions of the 
Glauber-Sitenko theory and tells us that the correspond- 
ing pd are no longer reliable. We use also a Gaussian 
ansatz for representing the amplitudes generated by the 
NN models A and D in analytical form. However, we 
aim at an excellent reproduction of the original ampli- 
tudes over the whole angular range. This means that 
typically we have to use 10 or even 12 terms in the sum 
in Eq. (TlOl) . instead of n = 5 which was taken in [25]. 

We examine the sensitivity by varying the large-angle 
behavior of the NN amplitudes that enter our Glauber 
calculation for pd. In a first scries of calculations a 
smooth cutoff on the angular region is introduced by 
multiplying the employed pN amplitude by the factor 
F(q) = 1/[1 + (? 2 A?o)1- Tms is done in the evaluation 
of the SS mechanism as well as for the DS mechanism, 



i.e. in the two-dimensional integral in Eq. ([S]). The lowest 
value for the cutoff momentum qq is chosen in such a way 
that it corresponds to a pN cm. scattering angle which 
is close to the position of the minimum of the elementary 
pN cross section. For the case of 179 MeV, shown in 
Fig. [3] (see dashed line), this amounts to 9 c . m . rs 105°. 
This choice for go makes sure that at least the pN cross 
section produced by those modified pN amplitudes shows 
a clear diffractive behaviour. We then increase the cutoff 
momentum qo until the corresponding angle 9 c .m. reached 
180 degrees. In the actual calculations the exponent v 
was varied in the interval v — 10 ™ 15. 

As an alternative to variations of go we performed also 
calculations with a sharp cutoff at the maximum trans- 
ferred momentum in the physical region q = 2kpN, where 
kpN is the pN cm. momentum. 

Finally, we employed different Gaussian parameteriza- 
tions, varying the number of terms between 6 and 12, 
and without cutoff. Those parameterizations of the pN 
amplitudes differ from the ones considered above in the 
angular range 90° < c . m . < 180° but, more importantly, 
for larger transferred momenta q outside of the physical 
range, i.e. q > 2kpN. Note that the integration in Eq. © 
requires the pN amplitudes at any (large) q, however, as 
already said above, the deuteron wave function strongly 
suppresses contributions to the integral from that region. 

The result of our analysis is summarized in Figs. 0H9] 
The bands represent the variation of the calculated pd ob- 
servables due to the cut-off procedures described above. 
We consider these bands as a sensible guideline for es- 
timating the angular region where the Glauber-Sitenko 
theory is able to provide solid results for a specific ob- 
servable and where this approach starts to fail. In partic- 
ular, they indicate when contributions from large angles 
start to become significant. Since such contributions are 
in contradiction with the basic approximations under- 
lying the Glauber-Sitenko model, any sizable influence 
from them undoubtedly marks the breakdown of this ap- 
proach. 

Based on our results we conclude that within the 
Glauber-Sitenko approach reliable predictions can be ob- 
tained for the differential cross section (Figs. 3HH1) an d 
also for the spin observables A y , A?, A xx , and A yy (Figs. 
[Jj [9]) for cm. scattering angles up to 50° — 60° in the 
pd system. Obviously, within this angular region there 
is practically no sensitivity to the pN amplitudes in the 
backward hemisphere, in accordance with the require- 
ments of the Glauber-Sitenko approach. As expected, 
for larger angles where such a sensitivity is observed, the 
width of the corresponding bands are smaller for higher 
energies (Fig. [5J and larger at lower energies (Fig. [6J. 
This feature can be seen in case of the differential cross 
sections and also for the spin observables A y , A yy , A xx 
presented in Figs. [7JIE1 and|Hl 

A comparison of our calculation with experimental 
data is only meaningful for such pd scattering angles 
where the Glauber-Sitenko approach works well. Ac- 
cording to our calculations this region includes the whole 
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FIG. 7: Spin observables of elastic pd scattering at 179 MeV versus the cm. scattering angle: Ay (a), A? (b), A yy (c), and 
Axx (d). Results of our full calculation (including the SS+DS mechanisms) are shown based on the NN models model A 
(green/grey) and D (red/black). The bands represent the sensitivy to variations of the large-angle tail of the pN amplitudes 
as discussed in the text. 



diffractive peak in the differential cross section da/dt at 
forward angles, for energies from around 50 MeV up- 
wards. This finding is important for the issue of the 
polarization buildup of antiprotons, because it validates 
the application of the optical theorem for evaluating the 
total polarized cross sections based on the obtained for- 
ward pd amplitude. 

With regard to the measured differential cross section 
at 179 MeV, see Fig. HI our Glauber-Sitenko calcula- 
tion describes quite well the first diffractive peak - for 
pN amplitudes generated from model A as well as for 
those of model D. The first minimum in the differen- 
tial cross section, located at q 2 m 0.12 — 0.13 (GeV/c) 2 
(i.e. 6 c . m . ~ 55°), and the onset of the second maximum 
is explained only by model D. The obvious strong dis- 
agreement with the data at larger transferred momenta, 
q 2 > 0.15 (GeV/c) 2 , corresponding to 8 c . m . > 60°, lies 
already in the region where the Glauber-Sitenko theory 
cannot be applicable anymore - for the case of the spin- 
less approximation for the pN amplitudes as well as when 
the spin-dependent amplitudes are included - and, there- 
fore, no conclusions can be drawn. Note that at lower 
energies, specifically at 50 MeV, the first minimum lies 
outside of the region where the Glauber-Sitenko approach 
can be trusted. In this context let us also mention that 
the results shown in Fig. [2] were obtained without any 
cutoff. 



Some remarks on the spin-dependent observables pre- 
sented in Figs. 13 and |9] The results obtained for 
the vector analyzing powers AP y and A d y indicate a strong 
model dependence. In contrast, the tensor analyzing 
powers A xx and A yy exhibit a very similar behaviour 
for both models A and D. We found that the spin- 
independent amplitudes dominate the latter observables 
and the inclusion of the spin-dependent amplitudes has 
only a minor influence. Thus, the results obtained here 
for A xx and A yy seem to be quite robust up to scattering 
angles of 60° — 70°. When the spin-dependent terms of 
the elementary pN amplitude (Bjy, Cat, C' n , Gn, Hn) 
are excluded, then the vector analyzing powers A? &ndA y 
vanish. 

At 50 MeV the uncertainties in the considered spin- 
dependent observables increase dramatically for angles 
around 65°, in accordance with the strong variations that 
one sees in the differential cross section (Fig. [6]), and, 
therefore, we do not show those quantities beyond 70°. 

We looked also at the influence of the D-wave com- 
ponent of the deuteron on the obtained results. In the 
differential cross section the contribution due to the D- 
wave is rather small in forward direction, but increases 
with increasing scattering angle. For example, at 179 
MeV the contribution by the D-wetve amounts to around 
30% of the absolute value in the region of the first mini- 
mum. The tensor analyzing powers A xx and A yy are con- 
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FIG. 8: Spin observables of elastic pd scattering at 282 MeV versus the cm. scattering angle. Same description of curves as in 

Fig. m 



siderably reduced (by one order of magnitude) when the 
D-wave is neglected. Actually, these observables practi- 
cally vanish if, in addition, the spin-dependent terms of 
the elementary pN amplitude are omitted. 



B. Total spin-dependent cross sections 

The total pd cross section is defined by @ 

otot = a + (7 1 PP-P d + a 2 (PP.k)(P d .k) + ( 73^, (11) 

where k is the unit vector in the direction of the an- 
tiproton beam, P p (P d ) is the polarization vector of the 
antiproton (deuteron), and P zz is the tensor polariza- 
tion of the deuteron (OZ||k). The total unpolarized 
cross section <7o and the spin-dependent cross sections <Ji 
(i = 1,2,3) are calculated using the generalized optical 
theorem as described in Ref. @, see also the Appendix. 
The results of our calculations are presented in Figs. ITUl 

eh na and ei 

Results for the total unpolarized pd cross section are 
displayed in Fig. [10] together with experimental informa- 
tion |33l - [37| . Obviously the unpolarized cross section is 
described rather well by both Jiilich models A and D [llj 
within the SS+DS mechanisms (dashed and solid lines, 
respectively), while it is overestimated by ~ 10 — 15% 
within the SS approximation, exemplified in Fig. UOI onlv 



for model D (dash-dotted line). A similar result was ob- 
tained in the spinless approximation in Ref. [§:]. Tak- 
ing into account the double-scatterig mechanism leads to 
more sizable changes in the results for the spin-dependent 
cross section, relevant for the spin-filtering mechanism, 
especially for 02 ■ One can see from Fig. [12] that this 
cross section is reduced by roughly a factor of two when 
the double-scattering mechanism is included. For the 
cross section <j\ this difference is in the order of 10-15%. 
Note that a decrease of the absolute values of <7o, o\ and 
02, due to shadowing effects, of comparable magnitude 
was reported in Ref. [l3| in a calculation based on the 
Nijmegen pN amplitudes. 

The cross sections 00 , <j\, and 02 are more or less com- 
pletely determined by the pN at forward angles and, 
thus, can be reliably calculated within the Glauber- 
Sitenko approach. Indeed, the uncertainty bands turned 
out to be very small and, therefere, we don't show them 
in the figures. The tensor polarized cross section 173 van- 
ishes in the SS approximation. At low energies 25-50 
MeV this cross section is about few mb in magnitude. 
Unlike the other cross sections discussed above 173 turned 
out to be fairly sensitive to the values of the pN ampli- 
tudes at large angles, i.e. to the variations considered 
in sect. Ill A. Thus, there is a significant uncertainty in 
the predictions based on the Glauber-Sitenko theory as 
indicated by the bands. With increasing energy the cross 
section 03 decreases and is only about 0.3-0.5 mb above 
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FIG. 9: Spin observables of elastic pd scattering at 50 MeV versus the cm. scattering angle. Same description of curves as in 

pig. m 



800 p 

700- 
600- 
500- 
^400- 
300- 
200- 
100- 




200 



Kalogeropoulos (1 980) 
Burrows (1970) 
Carroll (1974) 
Hamilton (1980) 
Bizzarri (1974) 




400 600 
P lab (MeV/c) 



800 



FIG. 10: Total unpolarized pd cross section versus the an- 
tiproton laboratory momentum. Results of our full calcula- 
tion (including the SS+DS mechanisms) are shown based on 
the NN models A (dashed line) and D (solid line). The re- 
sults obtained for the SS mechanism alone based on model 
D is indicated by the dash-dotted line. Data are taken from 
Refs. [H13l- 
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FIG. 11: Total pd cross section a\ versus the antiproton ki- 
netic energy in the laboratory system. Results are shown 
based on the NN models A (green/grey) and D (red/black). 
Calculations for the SS mechanism alone are indicated by dot- 
ted (A) and dash-dotted (D) lines while the full calculations 
(SS+DS mechanisms) are given by the dashed (A) and solid 
(D) lines. 



100 MeV. As expected, at higher energies the sensitivity 
to the pN amplitudes at large angles decreases too. 

With regard to the influence of the D-wave component 
of the deuteron wave function on the total cross sections 



we found that its contribution to o\ and 02 is less then 
1% for both considered NN models. The total unpolar- 
ized cross section ao decreases by ~ 5% if the D-wave 
is neglected, which is roughly comparable to the weight 
of the D-wave in the normalization of the deuteron wave 
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FIG. 12: Total pd cross section 02 versus the antiproton ki- 
netic energy in the laboratory system. Same description of 
curves as in Fig. 1111 
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FIG. 14: Dependence of the longitudinal (P\\, panel a) and 
transversal (P±, panel b) polarization on the beam energy. 
Same description of curves as in Fig. 1111 The acceptance 
angle is 20 mrad. 
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FIG. 13: Total pd cross section 03 versus the antiproton ki- 
netic energy in the laboratory system. Results of our full cal- 
culation (including the SS+DS mechanisms) are shown based 
on the NN models model A (green/grey) and D (red/black). 
The bands represent the sensitivy to variations of the large- 
angle tail of the pN amplitudes as discussed in the text. 



function. The cross section 173, which is non-zero only 
if the double-scattering mechanism is accounted for, is 
very sensitive to the -D-wave component. If the D-wave 
is neglected, then 03 changes significantly and, specifi- 
cally, remains positive over the whole considered energy 
range. 



C. Polarization efficiency 

Finally, let us discuss our results for the polariza- 
tion efficiency, the quantities relevant for the polarization 
buildup. According to the analysis of the kinetics of po- 
larization [j| HH , the polarization buildup is determined 
mainly by the ratio of the polarized total cross sections 
to the unpolarized one @. Let as define the unit vector 
C = Pt/Pt, where Pt = P d is the target polarization 



vector which enters Eq. (|11[) . The non-zero antiproton 
beam polarization vector P^, produced by the polariza- 
tion buildup, is collinear to the vector £ for any directions 
of Pt and can be calculated from consideration of the ki- 
netics of polarization. The general solution for the kinetic 
equation for pp scattering is given in Ref. |5|. Here we 
assume that this solution is valid for pd scattering also. 
Therefore, for the spin-filtering mechanism of the polar- 
ization buildup the polarization degree at the time t is 
given by @, [33 



Pp{t) = tanh 



where 



/ {do ± P T [*i + (C • k) 2 a 2 ] } 



(12) 



(13) 



Here n is the areal density of the target and / is the 
beam revolving frequency. One should note that the 
tensor cross section 03 from Eq. (fTTj) docs not con- 



tribute to Q^* 



Assuming the condition! Q_ — fl™ 



« 



(fP"*+fi° ut ), which was found in Refs. 0,111] for pp scat 
tering in rings at n — 10 14 cm~ 2 and / = 10 6 c _1 , one 
can simplify Eq. (|12[) . If one denotes the number of an- 
tiprotons in the beam at the time moment t as N(t), then 
the figure of merit is Pp(t)N(t). This value is maximal 
at the moment to = 2r, where r is the beam life time, 
which is determined by 00 the total cross section of the 
interaction of the antiprotons with the deuteron target 
as 



(14) 



nfa 



To estimate the efficiency of the polarization buildup 
mechanism it is instructive to calculate the polarization 
degree Pp at the time to [39] . In our definition for o\ and 
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tT2, which differ from that in Refs. 0,113; we nn d 

Pp(t ) = -2P T — , if C-k = , 
Co 

P p (t ) = -2P T ^±^, if |C • k| = 1 . (15) 

For evaluating the polarization degree Coulomb effects 
are taken into account via the procedure described in 
Ref. [||. Thus the quantities a, (i = 0, 1,2) in Eq. (15]) 
are actually the sum of the hadronic cross sections, of the 
Coulomb cross section and of the Coulomb nuclear inter- 
ference terms. In the concrete calculation an acceptance 
angle of 20 mrad is used. 

The polarization degree Pp(to) for £ • k = 1 (Pm) at 
Pt = P d = 1 is shown in Fig. [Til versus the beam energy 
together with the results for £ ■ k = (P±). One can see 
that, in general, the polarization efficiency increases with 
increasing energy. For longitudinal polarization maximal 
values of about 10-15% are predicted above 150 MeV, see 
Fig. [T4k). The transversal polarization degree is smaller 
than the longitudinal one for both models A and D. Of 
course, and as was already pointed out in our earlier 
works [1, H3|, there is a significant model dependence 
in the predictions for both polarization cases. 

Obviously, the inclusion of the DS mechanism leads 
to a decrease of the longitudinal P|| as well as of the 
transversal P± polarization efficiences by about 20-30% 
as compared to the SS mechanism alone [30]. Neverthe- 
less, for both considered models the magnitude of the 
spin-dependent cross sections is still comparable or even 
larger than those for pp 0, [3^, H(| ■ In this context let us 
also mention that the values for the polarization degree 
we obtained are somewhat smaller than those presented 
in [l3|, based on the amplitudes of the Nijmegen NN 
partial wave analysis [l4[ . 



spin-independent amplitude (with a pronounced diffrac- 
tive behaviour) so that the former do not spoil the ap- 
plicability of the Glauber-Sitenko theory. This means, 
in turn, that the approach can be used savely for the 
calculation of the total spin-dependent pd cross sections 
via the optical theorem. With regard to the considered 
pd differential cross sections and vector- and tensor an- 
alyzing powers (A^, A?, A yy: A xx ) our investigation in- 
dicates that reliable predictions can be obtained for cm. 
scattering angles up to 50° — 60° in the pd system. For 
179 MeV, where data on the differential cross section ex- 
ist, this range covers the first minimum and the onset 
of the second maximum. Here our results based on the 
pN amplitudes of the Julich model D turned out to agree 
nicely with the experiment while model A overestimates 
the measured pd cross section at the minimum. 

The total polarized pd cross sections Oi (i = 1,2,3), 
and specifically, the polarization efficiencies are quanti- 
ties that are of interest in the context of plans to es- 
tablish a polarized antiproton beam via the spin-filtering 
method as proposed by the PAX collaboration. Corre- 
sponding predictions presented in this work, exhibit a 
sizable model dependence, reflecting the uncertainties in 
the spin dependence of the elementary pp and pn interac- 
tions. Still, for both considered models we find that the 
magnitude of the spin-dependent cross sections is com- 
parable or even larger than those for pp. Thus, our re- 
sults suggest that deuteron targets can be used for the 
polarization buildup of antiprotons at beam energies of 
100-300 MeV with similar and possibly even higher effi- 
ciency than pp scattering. Nonetheless, only concrete ex- 
perimental data on the spin-dependent part of the cross 
sections of pp and pd scattering will allow one to confirm 
or disprove the feasibility of the spin filtering mechanism 
for the antiproton polarization buildup. 



IV. SUMMARY 

In the present work we analyzed the role of the spin 
dependence of the pN amplitude in elastic pd scattering 
for energies of 50-300 MeV of the incident antiproton 
on the basis of the Glauber-Sitenko theory. In the actual 
calculations we utilized element ary p N amplitudes gener- 
ated from the Julich NN model [ljj. The S- and D-wave 
components of the deuteron were included into the calcu- 
lation and the single- and double-scattering mechanisms 
were taken into account. 

Since some of the spin-dependent amplitudes exhibit 
a non-diffractive behaviour we performed various test 
calculations in order to pin down the angular range 
where the Glauber-Sitenko theory can be reliably ap- 
plied. Thereby, it turned out that this approach works 
rather well for the region of the forward peak, over the 
whole considered energy region. Obviously, for the con- 
sidered pN models those amplitudes with non-diffractive 
character are fairly small as compared to the dominant 
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Appendix: Invariant amplitudes Ai and spin 
observables for pd elastic scattering 

The scattering matrix Mpd in Eq. (TJ) can be expressed 
via the 12 invariant amplitudes Ai in the following way 
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Mp d = (Ai+A 2 cr-n) + (A3 + A4cr-n)(S-q) 2 
+ (A 5 + A 6 a- • n)(S • n) 2 + (A 7 <r • k)(S • k) 
+ (A 8 <T.q)[(S.q)(S.n) + (S-n)(S-q)] 
+ (A 9 + A 10 cr • n)(S ■ fi) + (A n tr ■ q)(S ■ q) 
+ (A 12 <r ■ k)[(S ■ k)(S ■ fi) + (S ■ fi)(S • k)] . 

(A.l) 

Here S = (er p + <x n )/2 is the total spin of the deuteron 
and the definition of the unit vectors k, q, and n is given 
right after Eq. Q. In the coordinate system as chosen 
in Ref. (25[, with the axes e x = q, e y = n, and e z = k 
the differential cross section da / dt and analyzing powers 
Ap, Ay, A xx , and A yy take the following forms |2o| 

da/dt = S 

= |Ai| 2 + \A 2 \ 2 + ^{z + Re{2A\(A 3 + A 5 ) 

+ 2A*(A i + A 6 ) + A* 3 A 5 +AlA 6 }}, (A.2) 

where Z-^l^l 2 . 

A% = 2Re[2(A* 1 +A* 3 +A* 5 )(A 2 + A i + A 6 ) 

+ A\A 2 - A* 3 A 6 - A* 4 A 5 + 2A* 9 A w ]^-\ (A.3) 

Ay = 2Re[(2Al+A* 3 + 2A* 5 )A 9 + (2A* 2 + Al 

+ 2A* 6 )A w +A* 7 A 12 + A;A 11 ]Z-\ (A.4) 



Ayy = {2(\A 5 \ 2 + \A 6 \ 2 + \A 9 \ 2 + \A 10 \ 2 ) 

- (|v4 3 | 2 + |^ 4 | 2 + l^rl 2 + l^sl 2 
+ \A n \ 2 + \A 12 \ 2 ) + 

+ 2Re[Al(2A 5 - A 3 ) + A* 2 (2A 6 -A 4 ) 

+ A%A h + AlA 6 ]}T,-\ (A.5) 

A xx = {2(|A 3 | 2 + |A 4 | 2 + |A 11 | 2 + |Ai 2 | 2 ) 

- (\A 5 \ 2 + \A 6 \ 2 + \A 7 \ 2 + \A 8 \ 2 
+ |A,| 2 + |A 10 | 2 ) + 

+ 2Re[A* 1 (2A 3 -A 5 ) + Al(2A 4l - A 6 ) 

+ A* 3 A 5 + AIAq}}^- 1 . (A.6) 



The amplitudes Ai can be rewritten as linear combi- 
nations of the amplitudes Fi used in our previous paper 
Q. The latter are defined within a different basis as 
compared to Eq. (|A. 1|) . In terms of the amplitudes Fi 
the observables given in Eqs. (|A.2[) - (|A.6I) coincide with 
those given in the Appendix of Ref. 8]. For collinear 
kinematics (q = 0) the scattering matrix Mpd(0) con- 
tains four independent terms 4 if and can be written as 

m 



Mp d . a p(Q) = giSap + (g 2 - gi)k a kp + 

igzVitapi + i{gi - g3)<Jtk i k j e al 3 J , (A. 7) 
where cr$ (i — x, y, z) are the Pauli spin matrices act- 
ing on the spin states of the antiproton, e a( g 7 is the 
fully antisymmetric tensor, k a are the Cartesian com- 
ponents of a unit vector k pointing along the beam 
momentum, and gi (i = 1,...,4) are complex ampli- 
tudes. When considering Eqs. (jA.ip and (|A.7[) together, 
one can find that in collinear kinematics with OZ||k: 
A 2 = A A = A 6 = A s = A 9 = A12 = 0, A 3 = A 5 , 
Aiq = An. The independent amplitudes gi are related 
with the amplitudes Ai via 



gi = Ai+ A 3 , g 2 = Ai+ A 3 + A 5 , 

g 3 = -A w , g 4 = -A 7 . (A.8) 



Taking into account the proper normalization of the scat- 
tering matrix Mpd and using those relations between the 
Ai's and gi's, the total hadronic cross sections in Eq. (jTTJ) 
given in Ref. [|[ on the basis of the generalized optical 
theorem can be rewritten as 



cr = 2spnlm{gi + g 2 ) , 

a i = A^Im(g 3 ) , 

cr 2 = 4y/7rlm(g 4 - g 3 ) , 

03 = k\/irlm(g2 - g 3 ) . (A. 9) 
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